Kohn-Luttinger superconductivity in graphene 
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We investigate the development of superconductivity in graphene when the Fermi level becomes 
close to one of the Van Hove singularities of the electron system. The origin of the pairing instability 
lies in the strong anisotropy of the e-e scattering at the Van Hove filling, which leads to a channel 
with attractive coupling when making the projection of the BCS vertex on the symmetry modes with 
nontrivial angular dependence along the Fermi line. We show that the scale of the superconducting 
instability may be pushed up to temperatures larger than 10 K, depending on the ability to tune 
the system to the proximity of the Van Hove singularity. 
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Since the fabrication in 2004 of single atomic layers 
of carbon, this new material (so-called graphene) has 
been attracting a lot of attentionfl]. The undoped sys- 
tem has conical valence and conduction bands meetingat 
two different Fermi points (known as Dirac points) [2, Q- 
This peculiar dispersion has shown to be at the ori- 
gin of a number of remarkable effects, like the exis- 
tence of a minimum conductivity at the charge neutrality 
point 0, Si, 01 . 

From the point of view of possible applications, the in- 
terest in graphene has been driven by the large electron 
mobilities attained in typical experimental samples. An- 
other remarkable property is that graphene can be used 
to build Josephson junctions when placed between super- 
conducting contacts [Sj. It becomes then quite intriguing 
whether graphene may support superconducting correla- 
tions on its own under suitable experimental conditions. 
On theoretical grounds, it is known that a model based 
on the conical dispersion requires a minimum strength of 
the pairing interaction for the development of a super- 
conducting instability in the undoped system [y|. There 
have been already several proposals to drive graphene 
towards a pairing instability upon doping, placing the 
emphasis on the role of topological defects [10|, the effect 
of a metallic coating 11], or the possibility of inducing 
superconductivity by electronic correlations [13, [l3|. 

In this paper we investigate a different route to super- 
conductivity in graphene, when the Fermi level is close to 
one of the Van Hove singularities (VHSs) of the electron 
system. These are points characterized by a divergent 
density of states, which has the effect of enhancing the 
magnetic and superconducting correlations 1J| . The ori- 
gin of the pairing instability lies in the strong anisotropy 
of the Fermi line at the Van Hove filling, and it pro- 
ceeds following in essence the same mechanism proposed 
by Kohn and Luttinger to show that superconductivity 
can arise out of purely repulsive interactions [l5[. This 
is possible as long as the e-e scattering becomes highly 
anisotropic, so that a channel with attractive coupling 
may appear when making the projection on the symme- 
try modes with nontrivial angular dependence over the 
Fermi surface. 

In the case of graphene, there are two VHSs located 
at about 3 eV above and below the Dirac points in the 
spectrum. Each of the singularities correspond to the 



existence of three inequivalent saddle points of the dis- 
persion at the boundary of the Brillouin zone, as shown 
in Fig. [T] In order to find the dominant electronic insta- 
bility arising from the divergent density of states, one has 
first to determine the shape of the Fermi line when the 
Fermi level is close to the VHS. For this purpose, we have 
characterized the energy contour lines around the saddle 
points of the valence band by means of a tight-binding 
model, suited to fit the dispersion e(k) known from angle- 
resolved photoemission spectroscopy (ARPES)[ig|. In 
the model, we have considered the transfer integrals for 
first, second, and third neighbors of the graphene lattice, 
labelled respectively by i, d and i', and the overlap inte- 
gral s between first neighbors. In terms of these param- 
eters, the Fermi velocity at the Dirac points is given by 
vp — (3/2) (t — 2t' + 3sd), which must be set to the value 
found in graphene, vp k, 2.7 eV. Moreover, the level of 
the saddle points relative to the Dirac points turns out 
to be 3d + (t — it' — 2d)/(l -I- s), which must correspond 
to the value « 2.7 eV found in ARPES [11]. With this 
input, we arrive at the two conditions 



t' w rf-2.7s 
t « 2.7 + 2d-5.4s-3sd 



(1) 
(2) 



Finally, we can adjust the parameters to reproduce the 
curvature of the dispersion at the saddle points [lg|, ar- 
riving at a dependence of the hopping d which is linear 
on s to very high accuracy: 



dw 0.07 + 2.8s + 0(s2 



(3) 



The ARPES dispersion around the saddle points can 
then be fitted leaving free the overlap integral s. The im- 
portant point is that the hopping parameter t' remains 
always constrained to a very small value, t' k, 0.1. This 
parameter controls the approximate nesting of the Fermi 
line, that is, the possibility of having regions in which 
the dispersion satisfies £(k) « — e(Q + k), with a fixed 
momentum Q. This is realized in the model with hop- 
ping parameters obtained from Eqs. ([1|)-([3|), as observed 
in a typical plot of energy contour lines shown in Fig. [T] 
The measure of the nesting is given by the susceptibility 
x(Q,a;) at the momentum Q connecting two inequiva- 
lent saddle points, which diverges in any event due to 
the singular density of states. We may approximate for 
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FIG. 1: Plot of energy contour lines around the saddle points 
of the valence band of graphene, obtained from a tight-binding 
model with first, second, and third-neighbor hopping param- 
eters given by Eqs. dU-® for s = 0.1 . 



instance the dispersion with the deviation (5k from the 
saddle point (27r/3,0) by 



e(k) « -a Ski + Ski 



Then we have 

x(Q,^) 



2\/37r2 a + (3 



log(cj) 



(4) 



(5) 



where the nesting instability at a = 3/3 appears in the 
dependence 



log 



1 + yZ/J/Sa 
1- ^y|3/3a 



log 



1 + yW/a 

1- vw^, 



(6) 



This has to be confronted with the susceptibility at van- 
ishing momentum transfer, which is 



x{o,^) 



1 1 

47r2 y^ 



log(w) 



(7) 



For the lower VHS in graphene, sensible values of the 
parameters are a « 7.41 and /3 « 2.03 (as obtained for 
s = 0.1), which give c' ~ 3.6. We see therefore that 
the susceptibility at momentum Q prevails over that at 
vanishing momentum. 

When the Fermi line becomes close to the saddle 
points, the large anisotropy of the e-e scattering may 
actually induce a pairing instability. This can be un- 
derstood within a renormalization group framework, by 
reelaborating the argument given by Kohn and Luttinger 
long time ago[15]. We will denote the interaction vertex 
by V{9, 9') for the particular case of BCS kinematics, in 
which the incoming (outgoing) particles collide with zero 



total spin and zero total momentum and the angle 6 [6') 
locates the position of the spin- up particle over the Fermi 
line. The BCS vertex gets corrections at low energies by 
the effect of the high-energy electron modes in slices be- 
tween energy A and A -I- dA about the Fermi level [rn|. 
The integration of these modes gives the variation 



dv{e,e') 



dK 

T 



dO" 9fc|| 1 
(27r)2 86" v{9") 



Vi9,9")Vi9",e') 



where 



(8) 
is the gradient of the dispersion and 
dk\\/d9" is the variation of the momentum along the slice 
parametrized by the angle 9" . We can write the above 
equation in more compact form by passing to the variable 



1 



^ d9" dh 



27m(A) Jo v{9") 89" 



(9) 



where the density of states n(A) is introduced so that the 
new variable also ranges from to 27r 18]. After defining 
the transformed vertex by V{(p, cj)') = V{9, 9'), we get 



dVi^,cj,') _n{A) 
(9 log A 27r 



2ir 



d<p"Viq^,r)V{<P",(b') (10) 



We can further decompose the vertex V{(j), (f>') in terms of 
the eigenmodes ^^^^(0) for the different representations 7 
of the point symmetry group. 



y((/., (/.') = 



E 



VZ,n'^7n{^)'i'Zi^') 



(11) 



We obtain then the set of coupled scaling equations 



a log A 



l{A) Y^VlsVs^n 



(12) 



In this framework, we recover the analogue of the Kohn- 
Luttingcr mechanism when any of the couplings V^^ „ 
turns out to be negative, in such a way that an unstable 
flow develops when the Fermi line is approached in the 
low-energy limit A — > 0. 

We remark that the scaling equation pO|) encodes the 
corrections to the BCS vertex that are logarithmically 
divergent at low energies in the particle-particle channel. 
If we start solving the scaling equation at an interme- 
diate energy scale A, the initial values of the function 
V{(j), (/)') will be dictated by the bare interaction as well 
as by regular corrections to it, given in general by finite 
diagrams in the particle-hole channel. Previous scaling 
analyses of electrons near a VHS have shown that only 
the momentum-independent component of the interac- 
tion potential is not irrelevant at low energies [l9J, which 
is consistent with the large screening effects from the di- 
vergent density of states. For this reason, we may con- 
sider that a bare on-site repulsion U between electron 
densities with opposite spin provides a sensible form of 
interaction close to the Van Hove filling. The relevant 





FIG. 2: Lowest-order particle-hole corrections to the BCS 
vertex for a bare on-site interaction. Full and dashed lines 
represent the propagation of electrons around different saddle 
points. 



point is that the bare short-range interaction gives rise 
to particle-hole corrections to the BCS vertex like those 
represented in Fig. [2l If we measure the angles cf), cf)' 
with respect to the x-axis, we observe for instance that 
the value of V(0,7r/3) will be enhanced by the particle- 
hole susceptibility (O at momentum Q, while y(0, 0) is 
enhanced by the particle-hole susceptibility ([7]) at zero 
momentum. The prevalence of x(Qj A) implies that 

V{0, 0) - y(0, ^/3) ^ 3C/2x(0, A) - 2C/2x(Q, A) < 

In general, we can anticipate the dominant terms in the 
modulation of the BCS vertex complying with the sym- 
metry of the Fermi line: 

V{(l), (/)) « Co + C6 cos(6(/)) + . . . (14) 

y(0, 0) « Co + C2 cos{2cj)) + c^ cos(40) + Cg cos(6(/)) + . .(.15) 
V{cj), TT - 0) w cj,' -t- c^' cos{2cj)) -f c" C0S(4(/.) 

+c[i' cos(6(/)) -f . . . (16) 

As we will see, several relations can be obtained between 
the coefficients in ([Til) - p^ by estimating the strength of 
the different scattering processes around the Fermi line. 
The Fourier expansions (fT4 |) - (fT6|) match well with the 
decomposition (fTTI) of the BCS vertex in terms of the ba- 
sis functions ^J^ ((/)). The point symmetry group is Cg^, 
which has 6 irreducible representations. Four of them 
are one-dimensional, with respective sets of basis func- 
tions given by {cos(6n0)}, {sin(6n(/))}, {cos((6ri -f 3)(/))}, 
{sin((6n -I- 3)^)} (n being always an integer). The 
other two representations are two-dimensional and have 
sets of basis functions which can be represented by 
{cos(r7T,(/)), sin(TO0)}, with the integer m running over all 
values that are not multiple of 3 and which are odd for 
one of the representations and even for the other. We can 
write therefore an expansion of the BCS vertex following 
(fTTjl and matching the modulations in P^ -(fTF 



1/(0, (/.') = K),o + 2^0,6 (cos(6(/.) + cos(6(/)')) 
-|-2V2,2(cos(20) cos(20') -f sin(20) sin(2(^')) 
-F2y2,4(cos(20) cos(40') - sin(2(^) sin(4<?!.') + (j) 
-^21/3,3 cos{3(j)) cos(30') + 2Vl 3 sin(30) sin(3(/)') 



One can check that the terms displayed in pT| ac- 
count for the dependence of the BCS vertex in Eqs. 
(|T4| -(fT6 |) . We identify actually the different coefficients 

2Vb,6, on the one 
= 4^0,6 + ^3^3, on 
can be translated 



c'2 - 21/2,2 + 2^2,4, c^ = 2F2,4,c[j = 

hand, and c^' = 4^2,4, c^' = 21/3,2, c^' 

the other hand. The constraint (fT3 

then to these couplings, since 1/(0, 0) — T/(0, 7r/3) is given 

by the combination 3(c2 + c^)/2 = 3(c2 + C4)/2. We find 

3(1^2,2 + 2^2,4) ~ 3C/2x(0, A) - 2C/2x(Q, A) (18) 

This already points at the existence of an unstable flow 
in the channel corresponding to the representation with 
d-wave symmetry as long as x(0. A) < x(Q, ^)- A closer 
inspection reveals actually that both couplings V2,2 and 
V2,4 must be negative. To show this, one more con- 
straint can be enforced by noticing that T/(7r/6, — 7r/6) — 
i/(7r/2, -7r/2) = 3(c^' - C4)/2. The relevant point is that 
the corrections to the BCS vertex at those angles are not 
singular at low energy A, as they involve scatterin g by 
particle- hole processes that scale at most as ~ -\/A ^0]. 
Therefore, the dominant contribution to both couplings 
V2,2 and V2,4 arises from the logarithmic dependence on 
A at the right-hand-side of Eq. (fT5|) . 

As long as V2,2 and 1/2,4 are negative, we have a pair- 
ing instability in the system, whose critical scale can be 
estimated by solving the coupled scaling equations ([12]) 
in the relevant symmetry channel. We can truncate the 
set of equations to 



d ( V2.2 
9 log A \V-^a 



V2A 
ViA 



n(A) 



1^2,2 

Va2 



V2A 
ViA 



^2,2 ^2,4 
^4,2 ViA 

(19) 
Eq. P^ can be easily solved by passing to the eigenval- 
ues Ai, A2 of the matrix of couplings. The scaling equa- 
tions read for them 



d 



a log A 



A, 



n(A)A' 



(20) 



We remark that V44 does not appear at the dominant 
level in the expansion of the BCS vertex, implying that 
such a coupling must be much smaller than those dis- 
played in (|17p . The two eigenvalues can be approximated 

then by Ai,2 



T/2 
''2,2 



4I/2M/2 



It is clear that 



(1^2,2 ± 

the positive eigenvalue Ai will vanish in the low-energy 
limit, while the signature of the pairing instability will 
be given by the growth of A2 towards very large negative 
values as A ^ 0. 

To find the behavior of the negative eigenvalue, we take 
for n{K) the density of states about a VHS at an energy 
/i from the Fermi level: 



n{e) 



1 



47r2 ^5^ 



log 



Ao 



Ml 



(21) 



.(17) 



In the above expression, Aq is an upper bound for the 
cutoff in our model of the VHS, that we take as 1 eV. 
The other important factor in the resolution of (pOI is the 
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FIG. 3: Plot of the temperature Tc obtained by identifying 
ksTc with the energy at which the eigenvalue A2 diverges, as 
a function of the deviation fj, of the VHS from the Fermi leveL 



choice of initial conditions for V2.2 and V2,4 at the upper 
cutoff Aq. In this respect, we have taken a value of C/ = 4 
eV for the bare on-site repulsion, which is between the es- 
timates made for graphite and carbon nanotubes[2l|. In 
order to go beyond the perturbative particle-hole correc- 
tions to the BCS vertex, we have summed up the series of 
leading logarithms obtained by iteration of the particle- 
hole susceptibilities in the diagrams of Fig. [2l Thus, we 
have constrained the initial couplings by the condition 



3(^2, 



2V2. 



U 



U 



1-3C/x(0,m) 1-2C/x(Q,m) 



(22) 



The precise values of V2,2 and V2,4 have been obtained 
by adding the other constraint mentioned below Eq. 
(PI), which reads 3(2^2,4 - ^2,2) = V^(7r/6, -7r/6) - 
V'(7r/2, — 7r/2). In these conditions, we have determined 
the point of the pairing instability in terms of the energy 
A at which the solution of Eq. (PD|) diverges, marking 
the position of a pole in the BCS vertex. The results are 
shown in Fig. [3l where we have represented the criti- 
cal point in a temperature scale after trading the energy 
variable A by the thermal energy ksT. 

We observe that the scale of the pairing instability de- 
pends drastically on the value of the chemical potential 



/x measuring the deviation of the VHS from the Fermi 
level. The plot of Fig. [3] shows anyhow that the in- 
stability exists irrespective of the value of /i. We recall 
in this respect that the Kohn-Luttinger mechanism was 
proposed to put forward the idea that any Fermi liquid 
is unstable at sufficiently low temperature [1 5|| . On the 
other hand, an inflection point can be seen in the plot 
of Fig. [3] for a value of fi slightly below 0.2 meV. That 
feature corresponds to the case in which the energy scale 
of the instability (and the scale of the gap) coincides with 
the deviation /i of the VHS from the Fermi level. Values 
of n to the left of the inflection point correspond there- 
fore to the regime where the pairing instability is driven 
all the way by the VHS. 

The possibility of finding a pairing instability in 
graphene at temperatures of the order of ~ 10 K may rely 
on the ability to make a fine tuning of the Fermi level to 
the VHS. This may be feasible as long as the proximity to 
the divergent density of states corresponds to a situation 
with a very large compressibility, which is energetically 
very favorable. It has been actually shown that the VHS 
may pin the Fermi level over a range of doping levels, 
forcing the breakdown of the uniform charge distribution 
into patches with uneven electron density (phase separa- 
tion) so that the Van Hove filling is reached in one of the 
phases [22i]. 

In conclusion, we have seen that placing graphene in 
the proximity of the VHS of its valence band may be 
a good instance to induce a superconducting instability 
in the electron system. The origin of this effect lies on 
the large anisotropy of the e-e scattering along the Fermi 
line, which leads to an attractive coupling in a channel 
with d-wave symmetry. We have shown that the scale of 
the pairing instability may be pushed up to temperatures 
larger than 10 K, depending on the ability to tune the 
system to the proximity of the VHS. 

The financial support of the Ministerio de Educacion 
y Ciencia (Spain) through grant FIS2005-05478-C02-02 
is gratefully acknowledged. 



[1] See, for instance, A. K. Geim and K. S. Novoselov, Na- 
ture Mater. 6, 183 (2007). [8] 

[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, [9] 

M. I. Katsnelson, 1. V. Grigorieva, S. V. Dubonos and A. 
A. Firsov, Nature 438, 197 (2005). [10] 

[3] Y. Zhang, Y.-W. Tan, H. L. Stormer and P. Kim, Nature 

438, 201 (2005). [11] 

[4] N. M. R. Peres, F. Guinea and A. H. Castro Neto, Phys. 

Rev. B 73, 125411 (2006). [12] 

[5] M. I. Katsnelson, Eur. Phys. J. B 51, 157 (2006). 

[6] J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz and C. [13] 
W. J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006). [14] 

[7] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 98, [15] 



076602 (2007). 

H. B. Heersche et al, Nature 446, 56 (2007). 

E. C. Marino and L. H. Nunes, Nucl. Phys. B 741, 404 

(2006). 

J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Phys. 

Rev. B 63, 134421 (2001). 

B. Uchoa and A. H. Castro Neto, Phys. Rev. Lett. 98, 
146801 (2007). 

A. M. Black- Schaffer and S. Doniach, Phys. Rev. B 75, 
134512 (2007). 

C. Honerkamp, Phys. Rev. Lett. 100, 146404 (2008). 
R. S. Markiewicz, J. Phys. Chem. Sol. 58, 1179 (1997). 
W. Kohn and J. M. Luttinger, Phys. Rev. Lett. 15, 524 



(1965). [20] J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Phys. 
[16] T. Ohta et ai, Phys. Rev. Lett. 98, 206802 (2007). Rev. Lett. 79, 3514 (1997) 

[17] R. Shankar, Rev. Mod. Phys. 66, 129 (1994). [21] E. Perfetto et al., Phys. Rev. B 76, 233408 (2007). 

[18] D. Zanchi and H. J. Schulz, Phys. Rev. B 54, 9509 (1996). [22] J. Gonzalez, Phys. Rev. B 63, 045114 (2001). 
[19] J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Nucl. 

Phys. B 485, 694 (1997). 



